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We propose an experimental protocol to study p-wave superfluidity in a spin-polarized cold Fermi gas tuned 
by an s-wave Feshbach resonance. A crucial ingredient is to add a quasi- ID optical lattice and tune the fillings of 
two spins to the s and p band, respectively. The pairing order parameter is confirmed to inherit p-wave symmetry 
in its center-of-mass motion. We find that it can further develop into a state of unexpected vr-phase modulation 
in a broad parameter regime. Measurable quantities are calculated, including time-of-flight distributions, radio¬ 
frequency spectra, and in situ phase-contrast imaging in an external trap. The yr-phase p-wave superfluid is 
reminiscent of the yr-state in superconductor-ferromagnet heterostructures hut differs in symmetry and origin. If 
observed, it would represent another example of p-wave pairing, first discovered in He-3 liquids. 


Coexistence of singlet s-wave superconductivity with fer¬ 
romagnetism is a long-standing issue in condensed matter 
physics iU]. One of the most interesting phenomena is the so- 
called TT-phase achieved in artificially fabricated heterostruc¬ 
tures of ferromagnetic and superconducting layers ffl, 
where the relative phase of the superconducting order param¬ 
eter between neighboring superconducting layers is tt. The 
TT-state offers new ways for studying the interplay between 
superconductivity and magnetism and has potential applica¬ 
tion for quantum computing in building up superconducting 
qubits through the tt phase shift HU. Different settings for 
its realization have been discussed, such as in high super- 
conductors inns and in spin-dependent optical lattices iflUl . 


In this letter, first we show that an unconventional p-wave 
TT-phase superfluid state emerges in the experimental system 
of a Fermi gas in a quasi-one dimensional optical lattice 11211 . 
Then we propose experimental protocols for observing this 
novel state by tuning the spin polarization. This is reminis¬ 
cent of the yr-state in superconductor-ferromagnet heterostruc¬ 
tures 111]. However, the tt phase shift of the superfluid gap here 
arises from a different mechanism—the relative inversion of 
the single particle band structures (s- and p-orbital bands) of 
the two spin components involved in the pairing. As a re¬ 
sult of this novel pairing mechanism, such a yr-phase super¬ 
fluid state has a distinctive feature—a center-of-mass (COM) 
p-wave symmetry, which distinguishes it from other yr-states 
in previous studies mu. We map out the phase diagram 
as a function of controllable experimental parameters—atom 
density and spin polarization. There is a large window for the 
predicted COM p-wave yr-phase superfluids in the phase dia¬ 
gram at low density and it occurs at higher critical temperature 
in relative scales, enhancing its potential for experimental re¬ 
alization. Striking signatures are calculated for experimental 
detection: (1) the shape of the density distribution in time- 
of-flight shows dramatic changes resulting from the pairing 
between different parity orbitals (i.e., s and p orbitals); (2) 
distinctive features are found in the density of states (DOS) 
such as the existence of the finite gap and midgap peak, which 
can be detected via radio frequency (rf) spectroscopy; (3) in 



FIG. I: Zero temperature phase diagram as a function of lattice filling 
and polarization when tp/ts — 8 , tp/ts = 0.05, t'g/ts = 0.05 
and U/ta = —9. pyrSF and pFFLO stand for different modulated 
COM p-wave superfluid states with the center-of-mass momentum 
of Cooper pairs located at Q = (yr/a, 0, 0) and Q 7 ^ {tt/a, 0, 0), 
respectively. NG-I refers to a normal gas (without pairing) where the 
|s t) band is fully filled while the |p F) band is partially filled. NG-II 
is another kind of normal state where the |s f) band is partially filled 
while the |p f) band is nearly empty. The grey area is forbidden due 
to the Fermi statistics constraint on the lattice filling. The grey line 
stands for the empty state. 


the presence of a shallow trap background, the density (in¬ 
cluding atom number density and pair density) distributions 
in real space exhibit shell structures composed of the pre¬ 
dicted p-wave superfluids. The orbital degrees of freedom 
play an essential role here; recently the research of higher or¬ 
bital bands in optical lattices has evolved rapidly 1 131] . For 
p-band fermions with attractive interaction, the chiral center- 
of-mass p-wave superfluidity in 2D il4|] and superfluids sim¬ 
ilar to the Fulde-Ferrell-Larkin-Ovchinnikov IlSll were found 
in theoretical studies. As we shall show with the model below, 
the pairing composed of different parity orbital fermions will 
lead to unexpected COM p-wave yr-phase superfluids. 


Effective model. Consider a Fermi gas with s-wave at¬ 
traction composed of two hyperfine states, to be referred to 
as spin f and f, loaded in a strongly anisotropic 3D cu- 
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bic optical lattice. In particular, we consider the lattice po- where M is an integer. Then, the Hamiltonian in the momen- 
tential Vql = X] sin^(fci,rQ) with lattice strengths turn space reads 


Vz = Vy ^ Vx, where fci, is the wave vector of the laser fields. 
In the deep lattice limit, the lattice potential at each site can 
be approximated by a harmonic oscillator and the lowest two 
energy levels are s and px orbital states for the lattice potential 
considered here, while other levels such as py and pz orbitals 
are well separated in energy. In the following, the Px orbital 
state is simply called the p orbital. Due to the strong confine¬ 
ment of the lattice potential in y and 2 directions, the system is 
dynamically separated into an array of quasi-one dimensional 
tubes. Further, let the gas be tuned with spin imbalance by the 
techniques developed in recent experiments 1 3- IS I. A key 
condition proposed here is to tune the imbalance sufficiently 
that the spin f and | Fermi levels reside in the s and p orbital 
bands respectively in order to hybridize the spin and orbital 
degrees of freedom. In the tight binding regime, the system is 
described by a multi-orbital Fermi Hubbard model 


H = -ts^Cl^{r)Cs^{r + ex) + tp^Cl^{r)Cpi{r + ex) 

r r 

- + ey) + C']^(r)C'si-(r -f e^)] 

r 

- tpJ2^C!l^^r)Cpiir + Cy) + C^^(r)C'p 4 ,(r -f e^)] -f h.c. 

r 

- Ft (r) - Fi 5Z 

r r 

+ ( 1 ) 


where ts and tp are the hopping amplitudes along the x di¬ 
rection for the s and p band fermions, respectively, while 
and tp are the hopping amplitude along the y and z directions. 
All the hopping amplitudes as introduced in Eq. O are posi¬ 
tive and the relative signs before them are fixed by the parity 
symmetry of the s and p orbital wave functions. C^rrij) is a 
fermionic annihilation operator for the spin cr component (f 
and 4 ,) fermion with the localized v (s and p) orbital located 
at the lattice site r, and is the corresponding chemical po¬ 
tential. The onsite interaction (last term in Eq. ([Til) is of the 
density-density type and arises from the interaction between 
two hyperfine states, which is highly tunable through the s- 
wave Eeshbach resonance in the ultracold atomic gases. Here 
we assume that the interaction strength is much smaller than 
the band gap. Therefore, the s-band fully filled spin down 
fermions are dynamically inert which are not included in the 
Hamiltonian (Eq. ([T]i). In this work, we focus on the case 
with attractive interaction where superfluidity is energetically 
favorable. 

Phase diagram at zero temperature. In order to study the 
superfluidity in our system, we apply the mean-field approxi¬ 
mation and assume the superfluid pairing is between different 
parity orbitals, i.e., between |s f) and \p T) states, in a general 

M 

form A(r) = 17(Cp^(r)C'st(r)) = ^ exp(zQm • r), 

m—1 


Hmf = 5]£«(k)Cl^(k)C,t(k) + ^£p(k)C^^(k)Cp;(k) 

k k 

M 

+ ^ ^(A;,Cpi(-k + Q^)at(k) + h.c.) 

m—1 k 

IA |2 

( 2 ) 
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where es(k) = —2tsCOskxa — 2t'g{coskya + coskzo) — 
p,-|-,ep(k) = 2tp cos kxa — 2tp{cos kya + coskza) — ^ 4 ,,k 
is the lattice momentum in the first Brillouin zone (BZ), N is 
the total number of lattice sites and a is the lattice constant. 

A fully self-consistent mean held calculation for the space- 
dependent order parameter is numerically challenging. We 
restrict our discussion to two forms of a variational ansatz. 
The first one is analogous to the Eulde-Eerrell state, where it 
assumes that M = 1 and the order parameter reads A(r) = 

A exp(iQ • r). The absolute value of superfluid gap is a con¬ 
stant, but the phase alternates from site to site. The second one 
is chosen by M = 2, Ai = A 2 = ^ ^nd Qi = —Q 2 , namely 
A(r) = Acos(Q • r), which is an analogue of the Larkin- 
Ovchinnikov state. This variational approach adopted here 
was previously justified by density-matrix-renormalization- 
group methods llQIl . Here we choose Q pointing along the 
X direction, say Q = Q(1,0, 0) to fully gap the Eermi sur¬ 
face of this quasi-one dimensional system. Eor EE like states, 
the mean-field Hamiltonian in Eq. (O can be diagonalized 
through the Bogoliubov transformation. The grand canonical 
potential of the system reads 

k 7 —it 

( 3 ) 

where C±(k) = 5 [es(k) - ep(-k-t-Q) ± 

•\/[es(k) -b ep(—k -f- Q)]^ -|-4|Ap] is the quasi-particle 
dispersion. While for the case of LO-like states, the Hamil¬ 
tonian can not be diagonalized analytically. However, the 
pairing terms link the one-particle states with momenta lying 
a single line in the BZ. Therefore, we can solve the resulting 
eigenproblem numerically for relatively large systems. 

Erom our numerics, we find that the free energy of the 
analogous LO states is always lower than that of the EE like 
phases, except at Q = (tt/o, 0, 0), where the EE and LO-like 
ansatz are equivalent. So the ground state of the system is 
a COM p-wave superfluid state with modulated pairing order 
parameter oc cos(Q • r), which breaks the translational sym¬ 
metry spontaneously. Qualitatively, that is because the ±Q 
pairing opens gaps on both sides of the Eermi surface, taking 
advantage of the available phase space for pairing, while the 
EE states only open a gap on one side. Since the dispersion 
of the p band is inverted with respect to that of the s band. 
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the pairing occurs between fermions with center-of-mass mo¬ 
mentum Q ~ kp-f + kpi, where kpf and kp^ are the two 
relevant Fermi momenta. When the occupation numbers of 
|s t) and \p 1) states are equal, the tt- phase superfuild state 
with Q = (tt/o, 0,0) is the ground state of the system. In real 
space, the pairing order parameter is a function of staggered 
signs along the x direction and obeys A(r) = — A(r -f ae^) 
when combined with the periodicity A(r) = A(r -f 2aex). 
The TT phase shift of the superfluid gap here arises from the rel¬ 
atively inverted single particle band structures directly, unlike 
in the conventional FFLO state ll^ . The predicted 7r-phase 
superfluid state is found to be quite robust. Even when the 
occupation number difference between | s t) and \p |) states 
is finite, the 7r-phase superfuild state is still the ground state. 
Particularly in the low density region rig + 1, there is a 

large window for this 7r-state. 

When the polarization p = is sufficiently large, 

the center-of-mass momentum will become incommensurate 
with the underlying lattice and this incommensurate COM p- 
wave state will be referred to as pFFLO. As shown in Fig.[T] a 
phase diagram as a function of atom density and polarization 
has been obtained. When the polarization is below a critical 
value, the system favors the 7r-phase superfluids. Beyond this 
value, a first order phase transition occurs from the 7r-phase 
to pFFLO superfluid states. It is worth to note that a large 
regime of parameters is found to exist in the phase diagram, 
making the experimental realization of the predicted new p- 
wave pairing phases simpler. 


Here we would like to stress two distinctions of the pre¬ 
dicted p-wave superfluid states with different types of modu¬ 
lated pairing from the conventional s-band FFLO 112 ill . One 
is that the Cooper pairs here are composed of different par¬ 
ity orbital fermions (i.e., s and p orbitals), which lead to the 
pairing order parameter with p-wave symmetry in the COM 
motion. The other is that the pairing order is modulated with 
the wavevector Q kp^ + kpi, rather than determined by the 
Fermi surface mismatch. 


Experimental signatures. The most distinctive feature of 
the predicted 7r-phase and pFFLO superfluids is that the pair¬ 
ing order parameters are spatially modulated and have p-wave 
symmetry in the COM motion. This leads to several charac¬ 
teristic experimental signatures. (A) the single particle mo¬ 
mentum distributions exhibit unique properties in the follow¬ 
ing two aspects. The first one is the shape of the density 
distribution in time-of-flight. We calculate the spin-resolved 
(or equivalently orbital-resolved) density distribution in the 
time-of-flight measurement assuming ballistic expansion as 

where k = mr/{ht), (j>vO^) is the Lourier transform of the 
i^-orbital band Wannier function and k = k mod G is 

the momentum in the first BZ corresponding to k (G is the 
primitive reciprocal lattice vector). As shown in Lig. |2] the 
highest peak for p band fermions is shifted from zero mo¬ 
mentum resulting from the non-trivial profile of the p-wave 
Wannier function superposed on the density distributions, dis- 


(a) (b) 



kxa/Ti: kx^/^ 


FIG. 2: Prediction of spin-resolved density distribution in time-of- 
flight, in (a) for |s f) fermions, while in (b) for \p F) fermions. 
The solid line shows the density defined in the main text along the 
kx-axis. The dash lines show the intensities of the Wannier orbital 
functions oc for comparison. Other parameters are = 

0.5, Upi = 0.5, tp/ts = 8, t'p/ts = 0.05, t'g/ts = 0.05 and U/ta = 
-9. 


Anguished from the s orbital fermions. The second aspect 
is a mirror-translational symmetry of the axial density dis¬ 
tributions of |s t) and |p F) fermions for the 7r-phase su¬ 


perfluid state. Lollowing the standard analysis 1112112211 . we 


f 


define the axial density distribution in the momentum space 
as nl^kx) = ^ j dkydk.{CU\x)C 

U<T (k)) for |s t) and 
Ip F) fermions, respectively. Lrom Eq. Q, the atom density 
distributions of |s f) and |p F) fermions can be calculated 
through nj/cr(k) = In the zero temperature, we And that 

n.,(k) = n,,(Q-k) = 

E(k) = [£:s(k)-Fep(-k -F Q)]^-F4|Ap and Q = (tt/o, 0,0). 
Therefore, the axial density distributions defined above satisfy 
the relation ng.f{kx) = n-p^itt/a — k^). It is also confirmed in 
our numerics as shown in Lig.[3|c). These signatures can be 
detected through polarization phase contrast imaging lll2ll . 


(B) the COM p-wave superfluid state here has the spatially 
nonuniform pairing order parameter, making it different from 
conventional p-wave superfluids. This leads to crucial dif¬ 
ference in Bogoliubov quasi-particle spectra. A finite en¬ 
ergy gap is shown in the spin-resolved (or equivalently the 
orbital-resolved) density of states (DOS) for the tt -phase su¬ 
perfluids (Lig. [^). Such spin-resolved DOS is calculated as 
P..{E) = i^J|<-|2^(L;-C„) + |ur|2b(£; + C„)], where 
{u'^^ v'^Y' the eigenvector corresponding to the eigenen- 
ergy C,n of the mean-field Hamiltonian Eq. (|2]) and the sum¬ 
mation runs over all the eigenenergy. This finite gap in the 
DOS gives direct evidence of superfluidity. Lor the pLLLO 
state, there is a midgap peak in the DOS, shown in Lig. [3b). 
This midgap peak results from the spatially modulated pair¬ 
ing order parameter which signifies the emergence of Andreev 
bound states lEsIl . The energy gap and midgap peak are found 
in the spin-resolved DOS for both |s f) \p i) fermions. 
Lor example, the DOS of |s f) fermions is shown in Lig. [3] 
Such spin-resolved DOS signatures can be detected via radio 
frequency (rf) spectroscopy ll23 - E^ . giving an experimental 
plausible probe of the predicted p-wave superfluids. 


(C) the predicted COM p-wave superfluids require neither 
spin-orbital coupling nor an induced second order effective 
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FIG. 4: Finite temperature phase diagram as a function of temper¬ 
ature and polarization at a certain lattice filling when tp/tg = 8, 
t'p/tg = 0.05, t'g/tg = 0.05 and U/tg = -9. 


FIG. 3: Top row plots density of states (DOS) pg^{E), (a) showing 
the finite energy gap for the vr-phase superfluid state and (b) showing 
the midgap peak for pFFLO resulting from Andreev bound states. 
Bottom row plots the axial density distributions of |s t) \p 4) 
fermions in momentum space for each of vr-phase (c) and pFFLO 
(d). The red and blue solid lines show n'^^{kx) and np^{kx) respec¬ 
tively, while the blue dots show np^{Q — k^). See main text for the 
definition of pg-^^E) and For the pFFLO state in (b), since 

there is a large polarization which can be considered as an effective 
external magnetic field, it leads to a shift of the density of states. 
Therefore, the midgap peak in (b) is not at S = 0. In (a) and (c), we 
choose Ug^ = 0.5 and npj. = 0.5, while in (b) and (d), ns-|- = 0.52 
and Upi = 0.45. Other parameters are the same as in Fig.[T] 



p-wave interaction. It arises directly from a purely s-wave 
two-body attraction. This leads to a significantly improved 
transition temperature, which is confirmed by our direct cal¬ 
culation of finite temperature phase transitions for the model 
Hamiltonian in Eq. ([T])— see Fig.|4]for the phase diagram. The 
transition temperature of the predicted COM p-wave super¬ 
fluids is indeed significantly improved as com pare d to other 
conventional relative p-wave superfluids 1127112811 . For in¬ 
stance, suppose a quasi-ID optical lattice has the lattice depth 
14 = 5i?r, Vy = Vz = ISEp and the wavelength of the lattice 
beams is 1064nm, the mean-field superfluid transition tem¬ 
perature can reach around 60nK when the s-wave scattering 
length between ®Li atoms is Os — 326ao 1^ . where ag de¬ 
notes the Bohr radius. Further increasing Ug, to around 600ao, 
the transition temperature can reach around 0.2/j,K, or even 
higher. 

Shell structure in a trap. In the following, we will discuss 
the effect of a harmonic trapping potential superposed on the 
optical lattices. Assuming that the harmonic trapping poten¬ 
tial is sufficiently shallow compared with the lattice depth, it is 
natural to apply the local density approximation (FDA) and let 
the chemical potential vary as a function of the position. Here 
we consider the trapping potential in the x direction. Fig. |5] 
shows various shell structures by our calculation. When the 
polarization is small, the region in the center of trap is the 
TT-phase superfluid, surrounded by a normal gas shell. While 
increasing the polarization, the region in the center of trap is 
no longer a superfluid, but a normal gas, surrounded by a pF- 


FIG. 5: Shell structures with a background trapping potential. The 
superfluid gap A and density profile n^t and npi of |s 4) and \p 4) 
fermions are shown as a function of the coordinate x, when tp/tg = 
8 , t'p/tg — 0.05, t'g/tg — 0.05 and U/tg = —9. The polarization 
P = is fixed at 0.3 and 0.03 for the first and second row, 

respectively. The frequency of the harmonic trap is chosen to be 120 
Hz. 


FLO superfluid shell, which in turn is surrounded by another 
shell of a normal gas. The density profile of |s f) and \p 4) 
fermions along the x direction are also demonstrated in Fig.|5] 
which can be detected through in situ phase-contrast imaging, 
providing an experimental plausible probe in a trapped atomic 
Fermi gas ll2|] . 

Let us take a specific experimental system, ®Li atoms, as 
an example. This Fermi gas is typically prepared with a to¬ 
tal atom number around 10^ in a trap 112111611 . In the recent 


experiment, it is loaded in an optical lattice and cooled down 
to the antiferromagnetic (Neel) transition scale lE^ . At the 
center of the trap, the density of one particle per site is main¬ 
tained. One can take advantage of this experimental devel¬ 
opment to achieve the 7r-phase p-wave state. One apparent 
method is to have a higher density, say, through tightening 
of the trap po tential or filling more atoms into the lattice ini¬ 
tially i30l432ll . Taking the peak density of atoms to be around 
2 X 10^^cm“^, the central filling of the lattice system is three 
particles per site when the lattice constant a = 532nm. An¬ 
other method is through tuning the polarization parameter of 
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the two hyperfine-spin components, so to populate one and 
only one of the two components of atoms into the p-orbital 
band. The two methods may be combined to optimize the ex¬ 
periment. 

Conclusion. We propose that the pairing between different 
parity orbital fermions can lead to a p-wave 7r-phase super¬ 
fluid state. The origin of the tt phase shift of the pairing order 
is distinct from the previous studies of 7r-states. We show that 
the predicted 7r-phase here occurs in a broad range in the phase 
diagram especially in the low density region. Increasing polar¬ 
ization, we And a phase transition from the 7r-phase state to an 
incommensurate COM p-wave superfluid. Signatures of the 
predicted p-wave superfluid states are calculated for time-of- 
flight and spectroscopic measurements in the regime of atom 
density and spin polarization deemed experimentally accessi¬ 
ble. The transition temperature and various shell structures 
with a background trapping potential are obtained for future 
experiments to explore these new forms of p-wave superfluid 
states. 
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